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ABSTRACT 

Using Suzaku observations of three neutron star low-mass X-ray binaries (Ser X— 1, 4U 1820—30 
and GX 349+2) we have found broad, asymmetric, relativistic Fe K emission lines in all three objects. 
These Fe K lines can be well fit by a model for lines from a relativistic accretion disk ('diskline'), 
allowing a measurement of the inner radius of the accretion disk, and hence an upper limit on the 
neutron star radius. These upper limits correspond to 14.5 — 16.5 km for a 1.4 M Q neutron star. The 
inner disk radii we measure with Fe K lines are in good agreement with the inner disk radii implied 
by kHz QPOs observed in both 4U 1820—30 and GX 349+2, supporting the inner disk origin for kHz 
QPOs. Additionally, the Fe K lines observed in these neutron stars are narrower than those in the 
black holes that are thought to be close to maximally spinning, as one would expect if inferences for 
spin arc robust. 

Subject headings: accretion, accretion disks — stars: neutron — X-rays: binaries 



1. INTRODUCTION 

The equation of state of the ultra-dense matter in neu- 
tron stars remains enigmatic, however, it is through accu- 
rate measurements of neutro n star radii and masses that 
it can be constr ained (e.g. lLattimer fe Prakasbl 120071: 
iKlahn et al.ll2006l ). Unfortunately, there are very few ac- 
cessible and robust observational measurements of neu- 
tron star radii. 

There are currently several methods that utilize X- 
ray spectroscopy/timing to determine neutron star radii. 
For instance, if one takes the origin of the kHz quasi- 
period oscillations (QPOs) to be associated with the 
inner edge of the accretion disk then this pla ces con- 
straints on mass and radius ([Miller et al.l |1998T >. How- 
ever, while these kHz QPOs can be measured pre- 
cisely, and are likely associated wi th the inner accretion 
disk, their origin i s not certain (jLamb fe Miller 1200 It 
Stella et al.l 11999b lAbramowicz et alJ 120031 : iTitarchukl 



20021) . Burst oscillations during thermonuclear X-ray 



burst s may be used to constra in the neutron star ra- 
dius (Bhattach aryya et al.ll200"5T ). but this method alone 
cannot provide very tight constraints as there are many 
unknown parameters. Alternatively, thermal continuum 
emission f rom the stellar surface can give a measure of ra- 
dius (e.g. Ivan Paradiis fe Lewinlll987l:lHeinke et alJl2006t 
IZavlinl 120071 : iHo et alJ l2007t iBogdanov et alJ 120071 ) but 
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suffer from distance and/or model uncertainties. Addi- 
tionally, X-ray absorption lines from the stellar surface 
would measure the gravitational redshift, but these lines 
appear to be extremely rare and difficult to observe, and 
the only detection i s of modest stat istical significance 
(jCottam et al.ll2002t iKong et al.ll2007fh In this paper we 
present a method for placing upper limits on n eutron 
star r adii using relativistic iron emission lines (sec Miller 
120071 for a review) that originate in the inner accretion 
disks around the neutron stars in X-ray binaries. 

Velocity shifts encoded into disk emission lines di- 
rectly reflect the radius at which the line is formed. In 
the X-ray spectra of accreting black holes and neutron 
stars, iron K-shell e mission lines a r e widely observed 
(iNandra et all 119971: iRevnoldd Il997t IWhite et all I1985L 
ll986tlHirano et al.lfl987l : lAsai et alll2000f ). The extreme 
red-wing of some iron emission lines implies that they 
are formed in the inner accretion disks and are primarily 
shaped by dynamics including relativistic Doppler shifts, 
due to the high velocities in the disk, and gravitational 
redshifts (hence the lines are referred to as being 'rela- 
tivistic'). Therefore, iron K lines can serve as incis ive 
measures of the inner disk radius ([Fabian 

These relativistic, asymmetric Fe K emission lines are 
well-known and well-studied in the X-ray spectra of 
bot h supermassive black holes in Act ive Galactic Nu- 
clei ([Nandra et al.lll997tlRevnolds!l 19971) and stel lar-mass 
black holes in X-ray binaries ( Miller et al.l 12004] ). In the 
case of neutron stars, iron emission lines are weaker, how- 
ever, and until lately prior observations have not clearly 
revealed a relativistic line profile, a nd the lines could 
be well fit by just a Gaussian (e.g.. IWhite et al.llT98 
Asai et alJl2000t lOosterbroek et aiTl2001t lDi~Salvo et al 



20051) . Only very recently have observations of Serpens 
X-l with XMM-Newton shown a relativistic, asymmetric 
line in a neutr on star low-mass X-ray binary (LMX B) for 
the first time (Bhattac haryya fe Strohmaverl l2007) . con- 
firming the inner disk nature of the lines in neutron star 
LMXBs. We note that tho ugh the iron line in Serpens X- 
1 has been observed before (IWhite et al.l [l986: Asa i et aTl 
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200 0i: lOosterbroek et al.ll2001[) . the new observations by 
Bhattacharvva fe Strohmayerl (|2007h were the first time 
a clear asymmetry has been seen in the line profile. 

The high effective area of the Suzaku X-ray telescope 
(jMitsuda et al.ll2~007l ) in the Fe K band combined with its 
broadband energy coverage, means that it is optimized 
for studies of inner disks and relativistic regimes using 
iron lines. Here we present Suzaku observations of the 
Fe K lines in three neutron star LMXBs (Serpens X-l, 
4U 1820-30 and GX 349+2). We find them all to show a 
relativistic asymmetric profile. In section [2] we describe 
the data analysis and results, and in section [3] we discuss 
the use of these lines as probes of neutron star radii, as 
well as comparing the results with kHz QPOs seen in 
these sources and with iron lines in black holes. Finally, 
we summarize our findings in section 3J 

2. DATA ANALYSIS AND RESULTS 

The data were obtained with the Suzaku X-ray tele- 
scope, which has several detectors on-board. The X-ray 
Imaging Spectrometer (XIS) detectors provide spectra 
from 0.5 - 10 keV, whilst the Hard X-ray Detector (HXD) 
PIN camera provides spectra at higher energies from 12- 
40 keV (with current calibrations). For observations of 
Serpens X-l and 4U 1820-30, Suzaku was operated at 
the XIS nominal pointing position, with the XIS detec- 
tors operating in 1 /4 window mode and 1-sec burst clock 
mode to prevent pile-up. The observations of GX 349+2 
were obtained from the Suzaku data archive. Two obser- 
vations were performed of approximately 25 ksec each, 
with the telescope pointing at the HXD nominal posi- 
tion. The XIS was operated with 1/8 window mode and 
0.3-sec burst clock mode. As the observations are only 4 
days apart and the source count-rate and spectral shape 
did not vary significantly, we add the spectra from both 
observations together to improve the signal-to-noise. The 
observations are summarized in Table [TJ 

The standard data spectral products for the XIS were 
used for all observations. We also used the canned re- 
sponse matrices, as custom responses cannot currently be 
created with the burst mode. For the HXD/PIN data, we 
extracted the spectrum from the cleaned PIN event files 
following the standard analysis threads. The PIN non-X- 
ray background was extracted from the time-dependent, 
observation specific model provided by the instrument 
team, and was combined with the standard model for 
the cosmic X-ray background. 

We fit the spectra using XSPEC v. 11 (lArnaudlll99l ). 
The continuum fit is performed from 1 — 9 keV in the 
XIS and 12 - 25 keV in the HXD/PIN. We ignore below 
1 keV as this is the region where the accumulation of 
some carbon-rich material onto the CCDs has changed 
the effective are curve, additionally we ignore the 1.5—2.5 
keV region where there are large calibration uncertainties 
due to an instrumental feature and the 4 — 7.5 keV region 
where the broadened iron line profile is present. We only 
fit the PIN spectrum up to 25 keV as this is the energy 
at which the background dominates over the source flux 
for these objects. While the XIS consists of four separate 
detectors, we only fit XIS 2 and XIS 3 as their calibra- 
tions appear to be most robust. The XIS 2, XIS 3 and 
PIN are all fit simultaneously. The broadband energy 
coverage provided by the combination of both the XIS 
and PIN allows an accurate determination of the contin- 
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Fig. 1. — Suzaku spectrum of Serpens X-l. For clarity only the 
XIS2 (1 - 9 keV) and PIN (12 - 25 keV) data are plotted, and the 
XIS2 data have been re-binned. The black line shows the over- 
all model, the colors distinguish the separate model components. 
Green is the blackbody, red the multicolor disk blackbody, orange 
the power-law and blue the diskline. 

uum shape either side of the Fe K band. This is vital in 
robustly determining the line profile. 

It is well known that neutron star spe ctra can be fit b y 
a wide variety of different models (e.g. iLin et al.l 1 2007) . 
For simplicity, we fit a model with two thermal compo- 
nents (multicolo r disk blackbod y + blackbody) as per 
a recent study (|Lin et al.l [2007') . In that study as the 
sources vary in luminosity, the measured temperatures 
of both thermal components in this model are found to 
follow L cx T 4 , suggesting a strong physical motivation, 
where as in other models the measured temperatures of 
the thermal component do not follow L cx T 4 as the 
source luminosity increases. An additional power-law 
comp onent is also required (as is also found in lLin et al.l 
2007), and we note here that all 3 spectral components 
are required. The spectral index of the power-law is al- 
lowed to vary between the XIS 2 and XIS 3 as this is re- 
quired to remove some slight differences in slope between 
the calibration of the different detectors. The power-law 
component for the PIN is tied to the XIS 2 value, and 
this is the value quoted in the Table [2] The parameters 
from the continuum spectral fits for each object are given 
in Table [2] and the spectral fits to Ser X— 1 are shown in 
Fig. Q] (spectral fits to the other objects are similar). 
All uncertainties quoted throughout the paper are at the 
90% confidence level. 

Fitting the continuum reveals significant iron lines in 
all three objects, which are all detected at greater than 
the 7-cr level in all objects (which we determine from 
an f-test). We show the lines in Fig. [2] Importantly, 
the lines in these objects are all clearly revealed as hav- 
ing broad, asymmetric profiles. While we tried fitting 
a variety of different continuum models, the model de- 
scribed above is the only model that successfully fits both 
the XIS and PIN data. However, various other models, 
such as a disk blackbody + power-law or a blackbody 
+ power-law model give acceptable fits to the XIS data 
alone. It is important to note the iron line profile recov- 
ered from these different continuum models is robust and 
changes very little in shape. This demonstrates that the 
line is not due to continuum components crossing in the 
iron band. 
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TABLE 1 

Summary of Suzaku observations 



Object 


Observation date 


XIS 


Exposure time of spectrum (ksec) 
PIN 


Serpens X— 1 


24 October 2006 


30.5 


27.5 


4U 1820-30 


14 September 2006 


17.9 


28.4 


GX 349+2 


14 March 2006 


25.3 


17.3 




19 March 2006 


28.1 


22.9 



TABLE 2 

Parameters of the continuum models 



Parameter 


Ser X-l 


4U 1820-30 


GX 349±2 


N H (10 22 cm' 2 


0.56 ±0.01 


0.26 ±0.01 


0.93 ±0.01 


Blackbody (keV) 


2.28 ± 0.02 


2.46 ±0.03 


2.22 ±0.01 


B-body normalization (10 — 2 ) 


4.9 ± 0.2 


3.2 ±0.1 


9.7 ±0.1 


Disk blackbody (keV) 


1.21 ±0.01 


1.15 ±0.02 


1.29 ±0.01 


Disk b-body normalization 


103 ±5 


55 ±3 


130 ±5 


Power-law index 


3.6 ± 0.4 


2.1 ±0.1 


1.4 ±0.1 


Power-law norm. (10 -2 ) 


18 ±3 


8.17 ±0.02 


2.52 ±0.03 


Flux a (10 -9 erg cm -2 s" 1 ) 


5.9 ± 0.3 


3.6 ±0.2 


10.2 ±0.5 


*2 (d.o.f.) 


1.09 (1936) 


0.99 (1936) 


1.17 (1936) 



Flux is evaluated in the 0.5-10 keV range 
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Fig. 2. — Iron lines in the three neutron star LMXBs observed 
with Suzaku. The ratio of count rate to continuum model is shown 
versus energy. Black points are for data from the XIS 2 detector, 
and red points are from the XIS 3 detector. The solid line is the 
best-fitting diskline model. 



For slowly rotating neutron stars, the gravitational po- 
tential around a neutron star is expected to be very close 
to a simple Schwarzschild potential. In fact, the neutron 
star spin for 4U 1820-30 and GX 349+2 inferred from 
the difference between the freqencies of the kHz QPOs is 
less than 300 Hz, which wo uld imply that the spin pa- 
rameter j = cJ/GM 2 < 0.3 ([Miller et al.lll998fl . Frame- 
dragging corrections for such spins sho uld give errors of 
much less than 10% (|Miller et al.lll998fh We therefore fit 
the iron line profiles with a model of line emission from 
a relativistic accretion disk a ssuming the Schwar zschild 
metric (the 'diskline' model. iFabian et al.l [T9891 . The 
line profiles in all objects can be well fit by this model. 
When fitting the lines we restrict the line energy to be- 
tween 6.4 — 6.97 keV (the allowed range for different ion- 
ization states of Fe K emission). The inner radius, in- 
clination and disk emissivity are free parameters in the 
fit, whilst the outer disk radius is fixed at a large value 
(1000 R G ). 

The inner radius, Ri n , of the accretion disk (in units of 
gravitational radii, Re = GM/c 2 ) is important for shap- 
ing the line, and is a parameter of the line fit. We are 
therefore able to measure the inner radius of the accre- 
tion disk in the three objects. The full 'diskline' pa- 
rameters are given in Table [3[ The inner radius of the 
accretion disk places an upper limit on the neutron star 
radius. In general, we find that these inner disk radii are 
in the range 7 — 8 Rg, which is equivalent to 14.5 — 16.5 
km for a 1.4 M Q neutron star. 

We checked that the 'diskline' is the appropriate line 
model to use. Alternative models have been developed 
for emission lines formed in the space-time around a spin- 
ning black hole. We find that fits with the 'laor' model 
(jLaorlll99lT ) for a maximally spinning black hole give in- 
ner disk radii that are consistent with the diskline fits. 
We chose the diskline model since reasonable expecta- 
tions for neutron star parameters suggest a spin param- 
eter < 0.3. 

In addition to fitting with the diskline model, it is also 



4 



Cackett et al. 



TABLE 3 

DlSKLINE PARAMETERS AND KHz QPO DETAILS 



Parameter 


Serpens X— 1 


4U 1820-30 


GX 349+2 


Line energy (keV) 
R,„(GM/c 2 ) 
Emissivity index, (3 
Inclination (degrees) 
Normalization (10 -3 ) 
Equivalent width (eV) 
Xl (d.o.f.) 


fi aq+0-15 

7.7 ±0.5 
-4 8 +0 - 4 
26 ± 2 
6.6 ±0.6 
132 ± 12 
1.05 (3851) 


6.97_o.i8 

6 7+ 14 
"•'-0.7 

-4 3+ 1 ' 1 
— 2.1 

19+ 5 
1.8 ±0.4 
51 ± 11 
1.00 (3851) 


6.97_o.o2 
8.0 ±0.4 
-4.1 ±0.3 
23 ± 1 
7.6 ±0.6 
76 ± 6 
1.06 (3851) 


Ri n (km) from Fe K line 
Upper kHz QPO (Hz) 
Rin (km) from kHz QPO 


15.9 ± 1.0 


13.8±^ 
1060 ± 20 
16.1 ±0.2 


16.5 ±0.8 
978 ±9 
17.0 ±0.1 



Note. — The Ri n measur ements given in km are evaluated at M = 1.4 Mq. Ser X— 1 does not have a detected kHz QPO. kHz QPO 
measurements are taken from Zhang ct al. ( 1998aj) for 4U 1820—30, where we choose the value where the upper kHz Q PO frequency reache s 
a maximum and remains constant with increasing count rate. For GX 349±2 we use the kHz QPO measurement from Zhang ct al. (1998b). 



possible to put a lower limit on the gravitational redshift, 
z = GM/Rc 2 , from fitting the asymmetric line profiles 
with two Gaussians, one narrow to fit the peak, and one 
broad to fit the red wing. We compare the measured 
center energy of the broad Gaussian component with the 
line rest energy. If the disk extends all the way down to 
the neutron star surface, one could naively interpret the 
energy shift as the surface redshift. Doing this we get 
rcdshifts in the range 0.1—0.2, providing a consistency 
check with 'diskline' fits. 

We also note that there is expected to be a 'boundary 
layer' between the inner accretion disk and the neutron 
star. This material is expected to be optically thick, and 
not in a Keplerian orbit around the star. The blackbody 
component of our continuum model can be interpreted as 
this component. In alternative continuum models this is 
often modelled as a Comptonized component, though the 
high effective optical depth means that it essentially be- 
haves as a blackbody. In neutron star spectra there does 
appear to be two thermal compo nents that follo w T 4 as 
the source luminosity increases dLin et al.ll200l . one a 
Keplerian accretion disk and the other a small, optically- 
thick boundary layer. Moreover, the line profiles we ob- 
serve are well fit by the 'diskline' model, suggesting that 
such a boundary layer does not strongly affect the line 
profile. In fact, the inner disk radii that we measure may 
imply that the boundary layer is small, for reasonable 
neutron star radii. Additionally, the frequencies of kHz 
QPOs observed in neutron star systems are likely associ- 
ated with the orbital frequency of the inner disk, and the 
frequency of this QP O is similar to the expected Keple- 
rian frequency there (|Miller et al.lll998h . indicating that 
the inner disk is likely in Keplerian orbits. 

Finally, we note that to achieve asymmetric lines are 
observed here the iron line emission must originate in the 
accretion disk. Alternative models for line broadening, 
such as a Comptonizing coronac would require a very 
high optical dept h, making the Comptoni zing material 
effectively a disk ([Reynolds fe Wilms! l200dh . 

3. DISCUSSION 

Suzaku observations of Ser X— 1, 4U 1820—30 and 
GX 349+2 have revealed broad, relativistic Fe K emis- 
sion lines in each of the objects. By fitting with the 
'diskline' model, we were able to measure the inner ra- 
dius of the accretion disk, and hence an upper limit 



on the neutron star radius. Only in one case previ- 
ously h as a clearly asymmetric line profil e been ob- 
served (|Bhattacharvva fe Strohmaverl I2007D . and that 
was in Ser X— 1 with XMM-Newton. The iron line 
that we see in Ser X— 1 with Suzaku is very similar in 
sh ape to that seen independent l y with XMM-Newton 
by iBhattacharvva fe Strohmaverl (|2007f ) , and the inner 
disk radii are similar, though not quite consistent within 
1-cr. The observations presented here are more sensi- 
tive and have a broader continuum coverage, allowing 
tighter constraints on the inner disk radius. The fact 
that relativis tic content is now clear in the ind epen- 
dent work of Bhattach arvva fe Strohmaverl (|2007f ) and 
in the more sensitive Suzaku spectra presented here, sig- 
nals that these lines are robust and can now be used to 
constrain the nature of neutron stars. 

In Fig. [3] we show the constraints on radius and mass 
from these new Suzaku observations. They are com- 
pletely independent of any previous constraints and rely 
on different physics. For comparison we include the con- 
straints from kHz QPOs observed in X-ray binaries. We 
also show mass-radius curves from a variety of equations 
of state. The curves shown are not meant to be exhaus- 
tive, but are just shown to demonstrate a range in mass- 
radius expected from different equations of state. Whilst 
the radius upper limits from our current observations 
cannot rule out any equations of state, the potential of 
using these lines can be seen. Moreover, if masses for 
these neutron stars are determined, then the region of 
allowed mass-radius space is significantly reduced. 

Indeed, our analysis of broad iron lines in two sources 
that exhibit kHz QPOs (4U 1820-30, GX 349+2; 
IZhang et al.lll998af bl) opens up a new way to estimate 
the masses of the neutron stars in these systems. The 
upper kHz QPO frequency is thought to be close to 
the orbital frequenc y at the inner edg e of the optically 
thick emission (e.g. iMiller et al.l I1998T ) . and thus from 
a radius close to the Ri n inferred from the iron line. 
In Schwarzschild coordinates the expressions for orbital 
speed and frequency measured at infinity are the same 
as their Newtonian counterparts, v m b = \J GM/r and 
^orb = (l/27r)-\/GAf/r 3 . Therefore, a simultaneous mea- 
surement of speed and frequency allows a measurement 
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Fig. 3. — Neutron star gravitational mass and radius constraints. 
The gray region is ruled our because of causality (the sound speed 
must be less than the speed of light). The orange region is the 
allowed region from the highest frequency QP O observed (1330 
kHz in 4U 0614+091, van Straaten et al. 2000) and assuming the 
neutron star spin, cJ/GM 2 = 0.2, where J is the stellar angu- 
lar momentum. The blue region is allowed from this work on 
neutron star iron lines - here we plot the most well-constrained 
upper limit, which is from GX 349+2. For comparison, we also 
plot the mass-radius curves for a variety of equations of state. 
All curves are as labelled in ILattimer &: P rakash ( 200l|) ( see re f- 
erences therein for details of the equat ions of state), except for 
NL4 which is from [A kmal ct al. (1998) and Z271 which is from 
Horowitz & Pickarcwicz (2001). 

of the stellar mass regardless of the radius: 



Frame-dragging corrections to this expression for a space- 
time with angular momentum parameter j = cJ/GM 2 
are only of order - j/(r/M) 3 / 2 (e.g., iMarkovid I2000D , 
which is less than 2% for the expected j < 0.3. Simul- 
taneous spectral and timing observations of 4U 1820-30 
would be especially valuable due to the tentative infer- 
ence of a mass M > 2 M from a saturation of QPO 
frequency vs. count-rate and color (Zhang et al. 1998a; 
Bloser et al. 2000; see Mendez 2006 for a dissenting 
view) . 

Pending such a coordinated analysis, in Table 3 we 
present a simplified consistency test between the timing 
and spectral data in 4U 1820-30 and GX 349+2. For 
a given QPO frequency, is, the radius is bounded from 

1 /2 

above by i? max = [GM/(47r 2 ^ 2 )] , hence for a given 
source the largest observed frequency yields the smallest 
upper limit i? max - We list these limits in Table |31 as- 
suming a mass M — 1.4 M©. The general consistency 
between spectral and timing limits is encouraging. We 
also see that the best values for 4U 1820—30 are some- 
what offset, and find that the agreement would be im- 
proved for higher mass (e.g., M — 2M© works), but the 
current broad errors on the iron line and lack of simul- 
taneous timing measurements prevent us from drawing 
conclusions at this time. 



We note that we find low values (< 30°) for the sys- 
tem inclinations in all cases. This may be a selec- 
tion effect as the observed objects were ch osen based 
on st rong lines seen in previous ASCA data (|Asai et al.l 
2000), and the strongest lines are present in systems 
with the lowest inclinations. Ser X— 1 and GX 349+2 
do not have previously well defined inclinations, though 
Ser X— 1 is expected to be < 60° as no energy depen- 
dent X-ray dips or eclipses are seen (|Franket alJflflSl 
and GX 349+2 is also suggested to have a low in clina- 
tion (|Wachter fe Margodll996l ; lO'Neill et al.ll2001h . For 
4U 1820—30, an inclina tion of 35 - 50° was d etermined 
from UV observations (jAnderson et al.l fl997t ). and ob- 
servations of the X-ra y superburst find limits th at are 
consistent with this (|Ballantvne fe Strohmav er 2004), 
though also allow slightly lower inclinations. Differences 
in the inclination may be reconciled if the inner disk is at 
a different inclination to the binary system. This may be 
possible as radio jets are not always seen to align with the 
binary inclination. For example, the radio jet in the black 
hole X-ray binary GRO J1655-40 (Hicllmin g fe Rupenl 
|1995|) is inclined at a different angle that the inclina- 
tion binary system determined from opti cal lightcurves 
(|Orosz fe BailvdfT997t [Greene et al.ll200lh . 

The presence of relativistic iron lines in neutron stars 
as well as in black holes suggests that despite the phys- 
ical and observational differences between the two types 
of objects, the accretion disks must be fairly similar. Ex- 
tremely broad iron lines in some black hole sources are 
thought to be due to spinning black holes allowing the 
inner disk to extend further in (e.g. iMiller et al.l 12004 
iBrenneman fe Reynolds! I2006D . The iron lines we ob- 
serve in these neutron stars strengthen this argument as 
the neutron star iron lines all have inner disk radii much 
greater than these extreme cases in black holes. 

Finally, in observations of iron lines in black holes, a 
strong reflection compone nt, or 'hump', is gene rally ob- 
served above 10 keV (e.g. iMiniutti et al.ll2007[) . In our 
observations we do not find strong evidence for such a 
reflection component. Currently, reflection spectra have 
only been calculated for an ionizing spectrum that is 
power-law in nature. We find, though, that neutron star 
spectra turn over and drop off faster than a power-law 
and so any reflection hump would be less apparent. 

4. CONCLUSIONS 

We have found broad, relativistic Fe K emission lines 
in 3 neutron star LMXBs (Ser X-l, 4U 1820-30 and 
GX 349+2) using Suzaku. From 'diskline' fits to the lines, 
we were able to measure the inner radius of the accretion 
disk, and hence place an upper limit on the neutron star 
radius. However, present measurements cannot currently 
rule out any reasonable equations of state. 

We found that the inner disk radii we measure with 
iron lines are consistent with the inner disk radii implied 
by kHz QPOs in 4U 1820-30 and GX 349+2, supporting 
the inner disk nature of kHz QPOs. Additionally the iron 
lines observed in these neutron stars are narrower than 
those in the black holes that are thought to be spinning, 
as one would expect if the spin interpretation is correct. 

Recently, new observations are yielding impor- 
tant progress in understanding neutron stars, such 
as co nstraining the allowed properties of the crust 
(e.g. ICackett et all 120061: IWatts fe Strohmaveij 120061 : 
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IStrohmaver fc Brown 2002). But, determining the equa- 
tion of state of matter in the core still remains a difficult 
problem. As currently one single method will not con- 
clusively determine the equation of state, it is critical to 
have multiple, independent methods to measure masses 
and radii of neutron stars. The advantage of using Fe K 
emission lines as a probe of the neutron star radius is that 
they only require short observations to clearly reveal the 
relativistic lines, and does not require any knowledge of 
the distance to the object. Additionally, there are likely 
to be many neutron star low-mass X -ray binaries that 
have these lines (e.g. lAsai et al.ll2000t ). Thus, using rela- 



tivistic iron lines in neutron stars offers a promising, com- 
plementary method to constrain neutron star radii. Fu- 
ture proposed X-ray observatories such as Constellation- 
X and XEUS will be able to place extremely tight limits 
on even weak lines. 



We thank A. Steiner and J. Lattimer for very kindly 
providing the mass-radius curves for various equations of 
state. JMM acknowledges funding from NASA through 
the Suzaku GO program. 
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